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Abstract. In colliding- wind binaries, shocks accelerate a fraction of the elec- 
trons up to relativistic speeds. These electrons then emit synchrotron radiation 
at radio wavelengths. Whether or not we detect this radiation depends on the 
size of the free- free absorption region in the stellar winds of both components. 
One expects long-period binaries to be detectable, but not the short-period ones. 
It was therefore surprising to find that Cyg 0B2 No. 8A (P = 21.9 d) does show 
variability locked with orbital phase. 

To investigate this, we developed a model for the relativistic electron gen- 
eration (including cooling and advection) and the radiative transfer of the syn- 
chrotron emission through the stellar wind. Using this model, we show that the 
synchrotron emitting region in Cyg 0B2 No. 8A does extend far enough beyond 
the free-free absorption region to generate orbit-locked variability in the radio 
flux. 

This model can also be applied to other non-thermal emitters and will prove 
useful in interpreting observations from future surveys, such as COBRaS - the 
Cyg 0B2 Radio Survey 



1. Introduction 



For most hot stars, the radio radiation they emit is due to thermal free-free 
emission by the ionized material in the stellar wind. A number of stars, how- 
ever, also show non-thermal radio emission, characterized by a negative spectral 
indejfl- This non-thermal emissi on is belie ved to be due to electrons that are 
Fermi-accelerated around shocks ( Bell 19781 ). As these relativistic electrons spi- 
ral around in the magnetic field, th ey emit synchrotroii radiation, which we de- 
tect as non-thermal radio emission (|Bieging et al.lll989l 1. In recent years, it has 
become clear that the shocks responsible for the Fermi-accelera tion in massive 



stars are those created in the colliding-wind region of a binary (jVan Loo et al 
20061 . Van Loo, these proceedings). 



Non-thermal radio emission from massive stars in general has been reviewed 
by Benaglia (these proceedings). In the present paper we will be mainly con- 
cerned with one speci fic colliding-wind binary, Cyg 0B2 No. 8A. From spec- 
troscopic observations, iDe Becker et al.l (2004, 20061 ) found it to be an 06If -|- 
05.5111(f)) system with a period of 21.908 ± 0.040 d and an eccentricity of 0.24 
=b 0.04. Other parameters relevant for our work are listed in Table [TJ 



^The spectral index is given by a, the exponent in the power law relation between flux and 
frequency: F^ (x v" . A stellar wind with a thermal spectrum has a ~ +0.6. 
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Table 1. Parameters of Cyg 0B2 No. 8A used in this work. Data frora 
iDe Becker etall (|2006D . 



primary secondary 



TeS (K) 


36800 


39200 


i?* (Rq) 


20.0 


14.8 


(Mq) 


44.1 


37.4 


log Lboi/^o 


5.82 


5.67 


M (M0yr-i) 


4.8 X 10-^ 


3.0 X 10-6 


Voo (km s"^) 


1873 


2107 



This binary has also been known to be a non-therm al radio emitt e r sinc e 
the first major survey of 0-type star radio emission by iBieging et al. ( 19891 ). 
It was recognized as such by its clearly negative spectral index. The colliding- 
wind nature of this object is further confir r ned b y information from its X-ray 
emission. Observations bv iDe Becker et al.l ( 20061 ) show an essentially thermal 
X-ray spectrum, but with an overluminosity of a factor 19-28 compared to the 
canonical Lx/Ltoo\ = 10"'' ratio for 0-type stars. This overluminosity is due to 
additional X-ray emission by material heated in the colliding-wind region. The 
X-ray light curve also shows variability at the level of ~ 20 %. Although the 
light curve is not well sampled, it suggests that the variability is phase-locked. 

Colliding-wind systems such as Cyg 0B2 No. 8A are important objects be- 
cause they can provide information on the Fermi acceleration mechanism, which 
is also relevant to other fields of astrophysics, such as interplanetary shocks and 
supernova remnants. The colliding-wind systems can also be used to constrain 
the amount of clumping and porosity in stellar winds. This is important for 
mass-loss rate determinations i n single stars, which are subject to considerable 
uncertainties ( Puis et al.l l2008l ) . Finally, the non-thermal radio emitters can 
help us detect binaries which are difficult to find spectroscopically and they 
are therefore useful for improving the determination of the binary frequency in 
clusters. 



2. Radio Variability 

Because of the importance of the colliding-wind binary Cyg 0B2 No. 8A, we 
decided to monitor its radio flux. We obtained a number of 6 cm continuum 
observations with the NRAC0 Very Large Array (VLA) during 2005 February 04 
to March 12, covering about 1.6 orbital periods. Figure [1] plots these data as a 
function of phase in the orbital period. The radio fluxes show clear variability 
and the variability is phase-locked with the orbital period. Furthermore, the 
data repeat well from one orbit to the next. Since these data were collected. 



■^The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 
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we have also reduced VLA archive data (Blomme et al. 2009, in preparation). 
These confirm the phase-locked variability and the orbit to orbit repeatability 
seen in Fig. [TJ 



o 

CD 




0.5 - 



0.0 



0.0 0.2 0.4 0.6 

phase (21.9-d period) 



0.8 



1.0 



Figure 1. 6 cm radio fluxes of Cyg 0B2 No. 8A, as a function of orbital 
phase. The filled circles indicates data from a first orbit, the filled triangles 
data from the subsequent orbit. 

It is, however, quite surprising that phase-locked variability is found in such 
a short-period binary. The main part of the colliding- wind region is expected to 
be comparable in dimension to the separation between the two stars. This size 
(of order a few stellar radii) is quite small compared to the size of the free-free 
absorption region of each of the components (of order 100 stellar radii). We 
therefore expected all the synchrotron emission from the colliding-wind region 
to be absorbed by the free-free absorption. Hence non-thermal radio emission 
fro m Cyg 0B2 No. 8A should have been undetectable. 



Blomme ( 20051 ) made a simple model for this, assuming the synchrotron 



emission to be a point source positioned between the two stellar components. 
Following the radiative transfer of the synchrotron emission through the free- 
free absorption region showed, as expected, no detectable level of non-thermal 
emission. As the data clearly contradict this, the effect of porosity was also in- 
vestigated. The porosity variant of the model assumes that the wind is clumped 
and that at least part of these clumps are optically thick. This allows radiation 
to escape more easily through the holes between the clumps, allowing us to look 
deeper into the stellar wind. Porosity models were indeed marginally capable 
of explaining the fact that we do detect non-thermal radiation from Cyg 0B2 
No. 8A. This conclusion is important, also in the context of single stars, as 
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it would help explain the important discrepancies fo und in the mass- loss rate 
determinations of massive stars (e.g. Puis et al. 20081 ). 



3. Modelling 

The point-source approximation for the synchrotron emitting region is cer- 
tain to be too simple. We here introduce a more sophisticated model for the 
colliding-wind region and its associated radio emission. The model does not 
solve the equations of hydrodynamics. Instead it defines the geometric position 
of the contact discontinuity between the two winds using the equations from 
Antokhin et al.l (j2004 ). On either side of the contact discontinuity a shock is 



formed. To simplify the model, we assume that these shocks and the contact 
discontinuity are at the same geometric place. 

At the shocks new relativistic electrons are generated due to the Fermi- 
acceleration mechanism. Their momenta follow a power-law distribution with a 
high-energy cut-off due to inverse Compton cooling. The fraction of electrons 
that become relativistic is determined by assuming that 5% of the shock energy 
goes into the electron acceleration. We then follow the electrons as they are 
advected outward along the contact discontinuity. We take into account that 
they cool down as they do so (due to inverse Compton and adiabatic cooling) 
and we follow them till they no longer emit synchrotron radiation at the radio 
wavelengths we are interested in. 

Next, we calculate the synchrotron emissivity along the contact disconti- 
nuity. We do this in a 2-dimensional model. As the situation is rotationally 
symmetric around the line connecting the two components, we can then rotate 
this 2-dimensional emissivity into a 3-dimensional simulation box. Free-free ab- 
sorption and emission are also added to the model. Finally, we determine the 
emergent int ensity and flux using the radiative transfer technique developed by 
Adam ( 1990l ). We calculate this for a number of phases along the orbit, thus 



determining the flux variability with orbital phase. The input parameters we 
use in this model are listed in Table [T l We also assumed an inclination angle of 
32°, as estimated bv lDe Becker et alT 



The top panel of Fig. [2] shows the resulting 6 cm radio fluxes. Clearly, this 
more sophisticated model is capable of providing variable non-thermal radiation, 
showing that at least part of the synchrotron emission can escape from the free- 
free absorption. The observed level of non-thermal radiation (about 1 mJy max) 
is not reached, but that could be due to one of the many assumptions in the 
model, such as the fraction of the shock energy that goes into the relativistic 
electrons or the magnitude of the magnetic field. 

The reason that this model does predict non-thermal radiation is that the 
synchrotron emission region is considerably larger than expected. This may 
seem surprising, but it is due to the fact that at larger distances from the 
central line connecting the two stars a much larger volume contributes to the 
synchrotron emission (as the situation is rotationally symmetric around that 
line). That partly compensates for the fact that the synchrotron emission per 
volume becomes weaker as we move away from the central line. The model also 
predicts that in high-resolution radio images where we can resolve the colliding- 
wind region we will still only see the central part of the synchrotron emission 
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region, as the emergent intensity of the outer part will be too low for it to be 
detected. In this way, th e model is consistent w ith the high-resolution radio 
observations (e.g. WR140. [Dougherty et al.ll2005l ). 
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Figure 2. Results from the more sophisticated synchrotron emission model. 
Top panel: the 6 cm radio fluxes predicted by the model, as a function of 
orbital phase. Bottom panel: the cumulative flux at phases 0.0 and 0.5, as a 
function of distance from the central line connecting the two stars. The largest 
contribution of the flux occurs between 1000-3000 Rq, i.e. approximately 10 
to 30 times the separation between the two stars. 

We next try to quantify the extent of the non-thermal emitting region, by 
determining where the major part of the detected non-thermal flux comes from. 
We do this by numerical experiment: we artificially set the synchrotron emission 
to zero beyond a certain radius and see how much non-thermal radiation is still 
detected. The bottom part of Fig. [2] shows that with a small cut-off radius we 
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only obtain about 0.12 mJy, which is the free-free emission of the two stellar 
winds. In the region 1000-3000 Rq we get the main contribution, with the 
flux levelling off beyond that. The major flux contribution therefore comes from 
distances of 1000-3000 Rq, i.e. approximately 10 to 30 times the separation 
between the two stars. 

Although the new model is quite successful in predicting non-thermal ra- 
diation, there is still one major problem: the phases at which maximum and 
minimum occur are incorrect. Observed maximum is around phase 0.1 and 
minimum at 0.7, but the predicted values are nearly exactly the reverse. 

The parameters listed in Table [J are based on ty pical values for the spec- 
tral types of both components ( De Becker et alll2006l ). Quite an error range is 



associated with these parameters, and we therefore explored if we could get a 
better fit to the observed data by changing some of the star and wind parame- 
ters. No exhaustive search of the parameter space was attempted, but in Fig. [3] 
we present an improved model, where the terminal velocity of the primary and 
the mass- loss rate of the secondary are increased, while the mass- loss rate of the 
primary is decreased. While the phase is still not correct, the difference with 
the observations is less. Furthermore, the general flux level is also higher and in 
better accordance with the observations. 
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Figure 3. Predicted 6 cm radio fluxes as a function of orbital phase for a 
model with changed stellar parameters compared to those listed in Table [1] 
Here we use Voo,i = 3500 kms-\ Mi = 2 x IQ-^ Mq yr'^ and M2 = 6 x 10"^ 
M0yr-i. 

To explain the remaining discre pancies, we explor e d a fu rther sophistication 
of the model. Following the work of lParkin Sz Pittard ( 200S), we iri troduced the 
effect of orbital motion. In the simplification bv lParkin &: Pittardl . this is done 
by assuming that the contact discontinuity consists of a "cap" in the inner 
part, which is formed instantaneously and a "ballistic" part in the outer part, 
where the material streams out further without any net force working on it. 
In the ballistic part, a "time-delay" effect needs to be taken into account: at 
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any given phase, the contact discontinuity in the balhstic part is the result of 
material emitted at the ends of the cap at some previous phase. Taking the 
orbital motion into account, this resu lts in a spiral-shaped fo rm of the contact 



aped lo i 

discontinuity in the orbital plane (see lParkin Sz Pittardf 200^ . their Fig. 10). 



Although the preliminary calculations presented at the meeting suggested 
that this would improve the agreement with the observations, we have so far 
not succeeding in finding a good model. As the synchrotron emission region in 
this model is wound up into a spiral, we detect a large part of it at any orbital 
phase and we therefore end up with considerably less variability as a function 
of orbital phase. We suspect that in the present case of Cyg 0B2 No. 8A, the 
con tact discontinuity is not as tightly wound up into a spiral as suggested by 
the Parkin Sz PittardI ( 20081 ) model. Some twisting of the contact discontinuity 



would undoubtedly be present, however, and this might be sufficient to explain 
the phase shift between the observations and simpler models. 

4. Other Stars 

We intend to apply the present model to a number of other binaries where 
we have a considerable amo unt of data on their non-t hermal radio emission, 
such as HD 168112 (Blomme et al.ll2005l 1. HD 167971 (Bl omme et "all [2007) and 



Cyg 0B2 No. 9 (jVan Loo et al.1 120081 ). We wih also extend the model to X-ray 



emission and optical spectra. 

We also expect a considerable number of colliding-wind binaries to be de- 
tected in future radio surveys. One upcoming survey is COBRaS, the Cyg 
0B2 Radio SurvejEI (PI: R-K. Prinja). This e-MERLIN legacy project has been 
awarded nearly 300 hrs to make a deep survey of the core of the Cyg 0B2 clus- 
ter. This cluster has already provided us with three examples of colliding-wind 
binaries (No. 5, No. 8A and No. 9). 

In our Galaxy, there are only a few examples of young, very massive clusters, 
of which Cyg 0B2 is one. The survey will make a mosaic of the core region at 
5 GHz (with a 1-sigma noise level of 3 //Jy) and at 1.6 GHz (1-sigma noise 
level of 7.5 //Jy). Its main science goals are to study (i) the mass-loss and 
evolution in massive stars; (ii) the formation, dynamics and content of massive 
OB associations; and (iii) the frequency of massive binaries and the incidence 
of non-thermal radiation The combination of the two frequencies will allow us 
to easily distinguish the non-thermal radio emitters from the thermal ones and 
thereby detect the colliding-wind binaries. The high sensitivity will result in the 
detection of a significant number of binary systems, allowing detailed statistical 
studies of the colliding-wind phenomenon. 



5. Conclusions 

Contrary to expectations, the short-period binary Cyg 0B2 No. 8A does show 
radio variability that is phase-locked with its 21.9-day orbital period. New mod- 
els for the synchrotron emission from the colliding-wind region show that this 



• ^http : //www . homepages . ucl . ac . uk/ ^ucapdwi/cobras/ 1 
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region is indeed extended enough that not all the emission is absorbed by the 
free-free absorption. At present, we still have the problem that the phases of 
minimum and maximum flux are not correctly modelled. We suspect that a 
combination of different wind parameters and some slight twist of the contact 
discontinuity due to orbital motion will be needed to explain these differences. 

We note that the present model does not require any porosity in the stellar 
winds of both components. Colliding-wind systems such as Cyg 0B2 No. 8A 
can therefore be highly relevant to interpreting single stars, where clumping 
and porosity have been proposed to explain the discrepancies in mass-loss rate 
determinations . 

In future, the present models will be applied to other non-thermal emitters. 
The model will also be extended to cover X-ray and optical wavelengths. It 
is expected that future surveys, such as COBRaS, will turn up a considerable 
number of non-thermally emitting colliding-wind binaries. 
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